Self-masking of spontaneous symmetry breaking in layer materials 



o 
o 

Oh! 
<■ 

O 

(N 



CO 



I 

C 

O 

o 



> 

(N 
(N 
(N 

VO 
O 
oo 
o 



Hiroyuki Yamase 

Max- Planck- Institute for Solid State Research, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 

(Dated: April 20, 2009) 

We study d-wave Fermi surface deformations (dFSD), the so-called Pomeranchuk instability, on 
bilayer and infinite-layer square lattices. Since the order parameter of the dFSD has Ising symmetry, 
there are two stacking patterns along the c axis: (+, +) and (+, — ). We find that, as long as the c axis 
dispersion is finite at the saddle points of the in-plane band dispersion, the (+, — ) stacking is usually 
favored independently on the details of interlayer coupling, yielding no macroscopic anisotropy. The 
dFSD provides unique spontaneous symmetry breaking that is self-masked in layer materials. 

PACS numbers: 71.18. +y, 71.10.Ay, 74.70.Pq, 74.72.-h 



Spontaneous symmetry breaking is a key concept in 
physics. Typical examples in condensed matter are crys- 
tallization, ferromagnetism, antiferromagnetism, ferro- 
electricity, and superconductivity. Recently a new type of 
symmetry breaking was found in the shape of the Fermi 
surface (FS) [Hi]. The FS usually fulfills the point-group 
symmetry of the underlying lattice structure. However, 
it was shown that such symmetry of the FS can be bro- 
ken by electron correlation effects in the two-dimensional 
(2D) t-J P, i, H and Hubbard 0, i, i, models on a 
square lattice: The FS expands along the k x direction 
and shrinks along the k y direction, or vice versa. This 
instability is characterized by a d-wave order parameter 
and we refer to it as d-wave Fermi surface deformations 
(dFSD) or d-wave Pomeranchuk instability. 

The dFSD state has the same symmetry as the elec- 
tronic nematic state Q. Implications of such reduced 
symmetry were widely discussed for high-T c cuprates @ . 
In particular, the strong a6-anisotropy of magnetic exci- 
tation spectrum observed for untwinned YBa2Cu30 y 
was well explained in terms of the dFSD and its com- 
petition with singlet pairing [TT1 |. Moreover many pecu- 
liar features observed for the bilayer ruthenate Sr3Rii2 07 
[l2l Il3j also turned out to be well captured in terms of 
the dFSD instability 0, M, H E3 • 

The dFSD is expected in layer materials, which are 
characterized by weak interlayer coupling. Since the or- 
der parameter of the dFSD has Ising symmetry, the in- 
terlayer coupling drives two possible stacking patterns 
between the layers, (+,+) and (+,—); we call the for- 
mer ferro-type (F) stacking and the latter antiferro-type 
(AF). In the latter case, macroscopic anisotropy does not 
appear, leading to self-masking of the underlying dFSD 
instability. 

In this Letter, we propose a microscopic theory for the 
dFSD phenomenon in realistic layer materials. We find 
that AF dFSD state is usually favored independently on 
the details of interlayer coupling. Thus the dFSD insta- 
bility provides unique spontaneous symmetry breaking 
that is self-masked macroscopically. This generic con- 
clusion is applicable to a wide range of layer materials, 
including recently discovered iron-based superconductors 

(HQ. 

We consider both bilayer and infinite-layer models. We 



first focus on a bilayer square lattice and analyze the 
following pure forward scattering model, 
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where c\^ a (c l k(7 ) creates (annihilates) an electron with mo- 
mentum k and spin a on the i — A and B planes; — 

So- c ko- c ko- 1S t ne number operator; N is the total number 
of sites on the i plane and /j, is the chemical potential. For 
hopping amplitudes t and t' between nearest and next- 
nearest neighbors on the square lattice, £k is given by 
£k = — 2t(cos k x + cos k y ) — At' cos k x cos k y . The forward 
scattering interaction, /kk' = —gd^d^, drives the dFSD 
instability, where dk = cos k x — cos k y and g > [13, HH ■ 
The last term in Hamiltonian ((T|) is the hybridization 
between A and B planes. We consider three types of c 



axis dispersions, ej; 
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the latter two dispersions vanish at 
k = (-7T, 0) and (0, 7r) while the former does not. This 
difference is crucial to the stacking of the dFSD. 

Hamiltonian ([l} is analyzed in the Hartree approxima- 
tion, which becomes exact in the thermodynamic limit 
in our model. We obtain the mean field 



k 



(2) 



which is nonzero only if the electronic state loses the four- 
fold symmetry of the square lattice and is thus the order 
parameter of the dFSD in the A(B) plane. The FS is 
elongated along the k x and k y directions for r] A ( B > > 
and r] A ( B ^ < 0, respectively, i.e., the order parame- 
ter has Ising symmetry. AF (F) stacking is defined by 
r\ A r\ B < 0(> 0). We determine the mean fields self- 
consistently under the constraint that each plane has the 
same electron density. We choose t'/t = 0.35 [TH, |22| 
and t z /t = 0.1 [lH, which is appropriate for the bi- 
layer ruthenate. £k has saddle points in k = (tt, 0) and 
(0,7r), and thus the van Hove energy of £k is given by 
fjP r, = 4t' = lAt. In the following we set t = 1 and shift 
the energy such that = 0. 
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FIG. 1: fp-T phase diagrams in the bilayer model with e k = —t z for several choices of g at t z — 0.1. Solid lines, T^ nd , denote 
second order transitions, while first order transitions are denoted by open circles, T c lst , and dotted lines, T lst ; the latter, present 
in panels (a) and (b), corresponds to a transition between F and AF; solid circles represent tricritical points. The insets magnify 
the regions around fi ~ —0.33 and T w 0.15 in (a), and around fi w 0.08 and T ~ 0.02 in (b). 



We first consider the c axis dispersion, eg = — t z , 
from which we can extract a generic conclusion about 
the stacking of the dFSD. Figure [TJa) shows the phase 
diagram obtained for g = 1. The rfFSD instability oc- 
curs around the van Hove energy of the in-plane disper- 
sion, i.e. /z=0 as in the case of the single-layer model 
f Ref . \2u) ; T c is almost unchanged by the presence of eg. 
The AF dFSD state, namely ry A = -g B ^ 0, is sta- 
ble in most of the region of the phase diagram. The F 
stacking appears in a very small region, as shown in the 
inset of Fig. [IJa). When the coupling constant g is re- 
duced, the F regions increase, close to the edges of the 
transition line. Yet the AF region is still much larger 
than the F ones, as it is shown in Fig. [Tf b) . Upon fur- 
ther reducing g, the F and AF regions split from one 
another [Fig. Q^c)] . After this "intermediate" phase, the 
property of the phase diagram drastically changes below 
g = 0.5, as shown in Fig. [TJd): No rfFSD instability ap- 
pears around /i=0, but the transition occurs around the 
van Hove energy of the bonding and antibonding bands, 
namely around /i v h = ±i z = ±0.1. The phase diagram 
eventually contains only F regions. This property does 
not change down to an infinitesimally small g. 

To gain an understanding of the dFSD stacking, we 
consider the Landau expansion of the grand canonical 
potential (per lattice site) with respect to r\ A and 77 s , 



J < 0, respectively, and is given by 
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where = £k — fi ± eg is the dispersion of the bonding 
and antibonding bands; /(A) = (e A ' T + 1) _1 is the Fermi 
function at temperature T and /' is its derivative. 

As long as a second order transition occurs around /i=0 
[Figs.QJa) and (b)], we can analyze Eq. ((U by expanding 
it with respect to eg and obtain 



where 



(5) 



(6) 



Since the k summation in Eq. © is dominated by the 
contribution from the saddle points of ek, we expand ek 
and eg around the saddle points, ek is then written as 
£k = t after choosing a suitable choice of momentum 
variables fc+ and fc_ with a cutoff momentum A, namely 
\k±\ < A; m is related to t and t'\ c?k is replaced by a 
constant, which just rescales eg, and we set = 1. 
For eg = —t z , Eq. © becomes 



The condition, X 2 — J 2 = 0, determines a second order 
transition. J drives AF and F stacking for J > and 
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FIG. 2: Coupling constant J between r) A and 77 s as a function 
of fj, at T = 0.2 (a) and 0.01 (b) for e£ = -t z . In (a), the 
function of —fi" 2 is also drawn. 



where e\ — At zero temperature Eq. d5|) reads J = 

— 2 ™* z . This /i~ 2 singularity is, however, cut off at a 

scale of the order of T°, the transition temperature at 
/i = 0. Therefore J becomes positive for < T°, as 
shown in Fig. [2ja). Because of the universal property 
of the dFSD instability obtained for a single layer model 
pll |. the dFSD phase is stabilized in the region < T c °. 
Therefore the sign of J is expected to be positive, leading 
to AF stacking of the dFSD in a major region of the phase 
diagram, as seen in Figs. Q] (a) and (b). 

If the coupling constant g is reduced, the temperature 
scale of the <iFSD becomes small and the energy scale 
of e k must be taken into account. We consider Eq. ((4]) 
instead of its lowest order expansion of Eq. ([5]). The 
scale of e£ is defined as the bilayer splitting at the saddle 
points of the in-plane dispersion, namely 



A z = |ej| at k = (it, 0) and (0, w) 



(8) 



in the present case and thus A z = t z . The positive J 
region is bounded by the scale A z , rather than T, as seen 
in Fig.^b). Then J becomes negative for |/z| > A z and 
shows steep dips at the van Hove energy of the bonding 
and antibonding bands, i.e. at /i v h = ±iz- This structure 
originates from the singular contribution of /' in Eq. Q. 
When the temperature scale becomes much smaller than 
A z , we have seen that the dFSD instability occurs around 
A 4 = MvH- Therefore we obtain F stacking, as seen in 

Fig. UK). 

So far we have considered a simple c axis dispersion. 
However, the above analysis holds also for other choices 
of a c axis dispersion, e.g. e£ = —t z (cosk x — cosfcj,) 2 , 
as long as A z ^ 0. Hence we obtain a generic conclu- 
sion for bilayer systems: F stacking prevails when the 
temperature scale of the dFSD is much smaller than A z , 
otherwise the major stacking is AF [23| . 

What happens for A z = 0? This can be seen from the 
analysis of Eq. (0. When e£ is expanded as e£ = (a\k + + 
b\k-)t z around the saddle points of the in-plane band dis- 
persion, the logarithmic singularity in Eq. ([7]) is weakened 
to e log |e| , leading to J « 8 ^J* gjli, The dispersion in 

this case may be taken as e£ = — 4i z cos % cos for 
which we obtain a\b\ < 0. The phase diagram indeed 
showed AF and F stacking for /1 < and /1 > 0. On the 
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FIG. 3: fi-T phase diagram in the infinite-layer model with 
= — t z cos ^- for t z = 0.1. The notation is the same as in 
Fig. 1. 
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other hand, when e£ is expanded as e£ = [a-2,k\_ 
c 2 k + k-)t z , we obtain J w 16 !"/ z (2a 2 b 2 + c|)log^-. 



Thus F (AF) stacking is stabilized for 2a 2 b 2 +c^ > 0(< 0); 
for e£ = — 2t z (cosk x + cosk y ), J becomes positive, indi- 
cating that the major stacking is F. As a conclusion, for 
A z = the stacking of the dFSD depends on the precise 
form of 

To gain more insight into the physics of layer materi- 
als, we now analyze an infinite-layer model. This is given 
by the same Hamiltonian {T]) except that t z is replaced 
by t z cos ^ and the k summation is taken also for k z ; 
here we assume two layers per unit cell. The conclusions 
obtained from Eq. (0) still hold, since the k z dependence 
appears only in e£ in Eq. ([6]) and just yields a factor of 
i after the k z summation. AF stacking is thus stabi- 
lized in a major region of the phase diagram for A z =/= 0. 
Unlike the situation in the bilayer model, bonding and 
antibonding bands now depend on k z . In this case, even 
if the temperature scale of the dFSD becomes very small, 
the instability still occurs around [i = [H| • The analy- 
sis of Eq. © then holds, yielding AF stacking [Fig. [5]. 

We have focused on a c axis dispersion, which is cer- 
tainly a relevant interlayer coupling in metallic systems. 
There may be other couplings such as Coulomb interac- 
tion £\ . V^nfnf and spin-spin interaction £\ . JySf ■ 

Sj*. However, these interactions do not yield a direct 
coupling between r\ A and T} B . They may just shift the 
chemical potential and renormalize the c axis dispersion, 
which is indeed the case in the Hartree-Fock approxima- 
tion, as long as the in-plane interaction /kk' is dominant. 

While our Hamiltonian (1) is simple, it is known that 
a forward scattering model works very well for Si'3Ru207 
[3 EH OS [13 and that such a forward scatteringinter- 
action turns out to be inherent in the t-J [H, H, 13] and 
Hubbard 0, [H, [|| 0] models, minimal models for cuprate 
superconductors. In addition, although the symmetry is 
different from the dFSD, a forward scattering model is 
also known to work for URu 2 Si2 [2a |. 

The bilayer ruthenate Sr-jRA^O? is a promising mate- 
rial for the dFSD instability [12|, LL3| • LDA calculations 
[22| show that the intrabilayer dispersion does not van- 
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ish at the saddle points k = (ir,0) and (0, 7r), that is 
A z ^ 0. Since the temperature scale of the dFSD is 
about 1 K and is much smaller than A z [23] , we expect F 
<iFSD stacking within each bilayer. On the other hand, 
there is a very weak interbilayer dispersion, which is ex- 
pected to have the form factor t' z cos cos ^ because 
the Ru sites shift by (i, |) between adjacent bilayers. 
The interbilayer dispersion thus vanishes at k = (tt, 0) 
and (0, tt), i.e. A 2 = 0. In this case, J is estimated as 
J ~ — as calculated above. This suggests AF and 
F stacking between adjacent bilayers, namely the stack- 
ing (++, ,++,--iJ--) and (++,++,++,++,•••) 

for /n < and fi > [2(|, respectively. 

High-temperature cuprate superconductors are also 
promising materials for the cZFSD instability [HI ■ While 
competition with singlet pairing may suppress the insta- 
bility, the ciFSD is still important, yielding strong cor- 
relations p], IH, which may lead to a giant anisotropic 
response to a small external anisotropy, e.g. lattice 
anisotropy. Hence the stacking is determined by the ex- 
ternal anisotropy and we expect F stacking for Y-based 
cuprates [llj and AF stacking in La-based cuprates in 
the low temperature tetragonal structure 

We have studied a generic tendency for the g?FSD 
stacking in bilayer and infinite-layer models. We have 
found that the stacking does not depend on the details 
of interlayer coupling as long as the c axis dispersion is fi- 
nite at the saddle points of the in-plane band dispersion, 



namely as long as A 2 ^ 0. For bilayer systems we have 
found F stacking only when the temperature scale of the 
dFSD is much smaller than A z [Fig.QJd)], and AF stack- 
ing otherwise [Figs, [lja)-(c)]. For infinite-layer systems 
instead the major stacking pattern is always AF [Fig.[3]. 

The dFSD instablity is a generic tendency in correlated 
electron systems as found in different theoretical models 
[UBS H3. The dFSD is thus expected for various ma- 
terials except if other instabilities prevail over the g?FSD 
instability. Since the condition A z 7^ is fulfilled for most 
materials unless there is a special symmetry reason, typ- 
ical materials with a high critical temperature for the 
dFSD instability may have AF stacking. While sponta- 
neous symmetry breaking is by definition recognized by 
symmetry breaking of the system, the dFSD provides a 
conceptually interesting property that spontaneous sym- 
metry breaking can be self-masked in realistic materials. 

Given the fact that antiferromagnetism was recognized 
about twenty centuries later than ferromagnetism (28l . 
I29I ] , it is not difficult to believe that the dFSD instability 
is hidden in various layer materials. Because of a coupling 
to the lattice, the easiest step to find AF g?FSD materials 
is to test a twofold lattice modulation along the c axis. 
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